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Abstract 

(N ■ 

An upper bound for the Wasserstein distance is provided in the general framework of the Wiener- 
Poisson space. Is obtained from this bound a second order Poincare-type inequality which is useful in 
terms of computations. For completeness sake, is made a survey of these results on the Wiener space, the 
Poisson space, and the Wiener-Poisson space, and showed several applications to central limit theorems 
\^ ' with relevant examples: linear functionals of Gaussian subordinated fields (where the subordinated field 

can be processes like fractional Brownian motion or the solution of the Ornstein-Uhlenbeck SDE driven 
by fractional Brownian motion), Poisson functionals in the first Poisson chaos restricted to "small" jumps 
(particularly fractional Levy processes) and the product of two Ornstein-Uhlenbeck processes (one in the 
Wiener space and the other in the Poisson space). Also, are obtained bounds for their rate of convergence 
to normality. 

1 Introduction 

In recent years many papers have looked at combining Stein's method with Malliavin calculus in order to 
£f) ' uncover new tools for proving varios central limit theorems (CLTs). For example, I. Nourdin and G. Peccati 
£> . derived an upper bound for the Wasserstein (Kantorovich) distance (and other distances) on the Wiener 
space using Stein's equation, |19j . Later, the same authors along with G. Reinert derived a second order 
Poincare(-type) inequality which is useful (in terms of computations) for proving CLTs, and which in fact, 
can be seen as a quantitative extension of the Stein's method from which upper bounds for the rate of 
convergence to normality can be found, [2U] . The first two authors with A. Reveillac extended these results 
to the multidimensional case, [21] . In [29], G. Peccati, J. L. Sole, M. S. Taqqu and F. Utzet, were able to find 
an upper bound, similar to the one in [19], for the Wasserstein distance in the Poisson space. G. Peccati and 
C. Zheng succeeded in extending this to the multi-dimensional case in [30j . All these works are important 
as they give quantitative tools for computing whether a random variable converges to normality or not, and 
if so, its rate of convergence. 



The upper bound inequality can be written as, 

d W (F,N) <E\l- (DF, -DL~ X F) U \ + E[(\x(DF) 2 \ , |DL _1 F|)«] 

where dw is the Wasserstein distance, N ~ A/(0, 1), D is the Malliavin derivative, L _1 is the inverse of the 
infinitesimal generator of the Ornstein-Uhlenbeck (0-U) semigroup, and % is the underlying Hilbert space. 
In the case of the Wiener space upper bound, D is the Malliavin derivative defined in that space so this 
inequality holds even when the underlying Hilbert space is not Lz. Also, since there are no jumps here, 
the second term on the right disappears (x is the size of the jump). In the Poisson space case, since we do 
not (as yet) have a Malliavin calculus theory developed for a general abstract Hilbert space, the underlying 
Hilbert space must be H = L\. A question naturally arises - can this be done for a general Levy process; 
that is, is this upper bound achievable in a mixed space, the Wiener-Poisson space? The main difficulty 
in answering this question is that in the Wiener-Poisson space we don't have decomposition in orthogonal 
polynomials (unlike with Hermite polynomials in the Wiener space, see |12] for a complete explanation), and 
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results like the equivalence between the Mehler semigroup and the Ornstein-Uhlenbeck semigroup. So, to 
overcome these shortcoming, will be necessary to deduce new formulas that will allow us to follow the ideas 
developed in jTO] and [3U], and recover their results for the Wiener-Poisson space. Will be shown that this 
bound still holds even when both spaces are involved. 

Before get into the details, some notation: The measure [i on the underlying Hilbert space L 2 is defined by 
the underlying Levy process, that is, let L t be a Levy process (L t has stationary and independent increments, 
is continuous in probability and Xq = 0, with -E[£ 2 ] < oo) with Levy-triplet given by (0, a 2 , v), where v is 
the Levy measure, then Jf R+xR f(z)d/i(z) = a 2 J R+ f(t, 0)dt + iJ R+xRo f(t, x)x 2 dtdv{x), where M = R- {0}. 
On the other hand, in order to define a Malliavin derivative in the Wiener-Poisson space it is sufficient to 
have a chaos decomposition of the space L 2 (f2). This was achieved in P3] by K. Ito, so any random variable 
in L 2 (£l) has a projection on the q th chaos given by I q (f q ), where f q is a symmetric function in L 2 . Also, 
when working in the Wiener-Poisson space, the Malliavin derivative can be regarded in terms of "directions" , 
i.e., we can think of it as the derivative in the Wiener direction or the derivative in the Poisson direction. 
The fact that this can be done in this way was proved in [SI] by J. L. Sole, F. Utzet and J. Vives (a quick 
review of the theory is given below). They showed that the Malliavin derivative with parameter z G R + x R 
can be split in two cases, when z = (t, 0) and when z(t, x) with x =/= 0. The first case will be the derivative 
in the Wiener direction (intuitively because there are no jumps when x = 0), and the second will be the 
derivative in the Poisson direction. A distinction between the Malliavin calculus in the Wiener space or the 
Poisson space and in this Wiener-Poisson space is the need to define two subspaces of L 2 (fl): one where 
the Malliavin derivative in the Wiener direction coincides with the usual Malliavin derivative in the Wiener 
space and is well defined, and another where the Malliavin derivative in the Poisson direction is well defined. 
These subspaces are denoted by DomD w and Dom£)' 7 . 



Theorem l.(Main result 1: upper bound) 

Let N ~ Af(0, 1) and let F eDomD w DDomD J be such that E[F] = 0. Then, 

d w (F,N)<E 1 — (DF, —DL~ 1 F) 2 +E[{\x{DF) 2 \ADL- 1 F\) 2 ] 

As mentioned above, in [20] , the authors used the respective Wiener space upper bound to deduce a second 
order Poincare inequality. In this sense, our second main result is the second order Poincare inequality in 
the Wiener-Poisson space derived from the above upper bound. This version is weaker than the one in the 
Wiener space because of the scarcity of results in this space, that would enable us to prove it for functionals 
in DomD M/ nDom£) J ; nevertheless, is possible to state it for functionals that lie in one specific Ito-chaos. 



Corollary l.(Main result 2: Second order Poincare-type inequality) 

Fix q £ N and let F = I q (f) with E[F] — /i and Var[F] — a 2 , and assume that N ~ A/"(/i, er 2 ). Then, 
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Thanks to these (Wiener, Poisson and Wiener-Poisson) upper bounds and second order Poincare inequal- 
ities, many CLTs can be proved and generalizations made. In this paper, these bounds are reviewed for each 
space, showing their importance by giving applications with relevant examples. In the Wiener space case 
the second order Poincare inequality is used to prove CLTs for linear functionals of Gaussian-subordinated 
fields when the decay rate of the covariance function of the underlying Gaussian process satisfies certain con- 
ditions. These CLTs are applied to the important cases where the underlying Gaussian process is either the 
fractional Brownian motion or the fractional-driven Ornstein-Uhlenbeck process, with H G (0,|) U (|,l). 

In the Poisson space case, the respective upper bound is used to prove that the small jumps process 
(jumps with length less than or equal to e) of a Poisson functional process with infinitely many jumps goes 
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to a normal random variable when e goes to zero. Furthermore, a remarkable extension, of the known result 
(proved in 1 ) which states that the small jumps process of a Levy process can be approximated by Brownian 
motion as e goes to zero, is proved. It is extended to Poisson functionals (ii(/)J and showed that the small 
jumps process of this functional can be approximated by a Gaussian functional with the same kernel / as 
e goes to zero. Then is applyed this result to show that in order to simulate a fractional (pure jump) Levy 
process (FLP), it is sufficient to simulate a process with finitely many jumps plus an independent fractional 
Brownian motion (FBM). 

Finally, the second order Poincare(-type) inequality, developed in this paper, is used to prove that the time 
average of the product of a Wiener Ornstein-Uhlenbeck process with a Poisson Ornstein-Uhlenbeck process 
converges to a normal random variable as time goes to infinity. This example highlights the importance of 
the inequality in the Wiener-Poisson space, since it cannot be achieved by the upper bounds in the Wiener 
or Poisson space. An estimate of the rate of convergence to normality is obtained in the examples where the 
second order Poincar inequalities is used. 

The paper is organized as follows: In Section 2 is recalled the basic tools of Malliavin calculus on the 
Wiener space as well as the second order Poincare inequality. Then is stated the Malliavin calculus results 
for the Wiener-Poisson space and explain how these tools have been extended. In Section 3 the theory 
developed in [20] and [19] is reproduced using Malliavin calculus theory for the Wiener-Poisson space, and 
with this framework is stated a "Levy version" of the second order Poincare inequality. Section 4 is dedicated 
to going over the inequalities for the Wiener, Poisson and Wiener-Poisson spaces. In the Wiener space case, 
a result proved in [3D] concerning CLTs of linear functionals of Gaussian-subordinated fields is extended. In 
the Poisson space case, a result on the simulation of small jumps for processes with infinitely many jumps 
is given. Finally, an example of application of the second order Poincare inequality in the Wiener-Poisson 
space is showed. 

2 Preliminaries 

As mentioned above, the most important tool for proving CLTs on the Wiener space is the upper bound 
developed in J5DJ. This requires various Malliavin calculus results on the Wiener space (Malliavin derivative, 
contraction of order r, Mehler formula, etc) which are extensively studied and explained in [25 . For the 
sake of the completeness, the basic tools from Malliavin calculus in the Wiener space are reviewed and then 
is defined the Malliavin calculus in the Wiener-Poisson setting, both needed in this paper. 

2.1 Malliavin Calculus on Wiener space 

Let H be a real separable Hilbert space. Assume there is a probability space (f^J 7 , P) where W := 
{W(h) I h £ H} is an isonormal Gaussian process, that is, W is a centered Gaussian family s.t. E[W (hi)W (/12)] = 
(hi,ti2)-H- Choose T to be the c-algebra generated by W. Let H q be the q th Hermite polynomial, 

H q {x) = {—l) q e~-^(e~~), and define the q Wiener chaos of W (denoted by H g ) as the subspace of 
L 2 (tt) := L 2 (Q,T,F) generated by {H q (W{h))/ h e H, \\h\\ n = 1}. Is important to underline that L 2 (VL) can 
be decomposed (Wiener chaos expansion) into an infinite orthogonal sum of the spaces H g : L 2 (f2) = ©^ H g 

Remark 1. In the case where H — L 2 ; for any F G L 2 (f2), 

F = £/,(/,) (1) 

q=0 

where I q is the q th multiple stochastic integral, fo — E[F], Iq is the identity mapping on constants and 
f q G are symmetric functions uniquely determined by F. 

Let S be the class of smooth random variables, i.e., if F E S then there exists a function (f> € C°°(K™) 
s.t. p^{x) has polynomial growth for all k = 0, 1, 2, 3, . . . and F — <f)(W(hi), . . . , W(h n )), hi € U. The 
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Malliavin derivative of F 6 S with respect to W is the element of L 2 (il, Ti) defined as 

DF = jr^-(W(h 1 ),...,W(h n ))h l (2) 

i—l 

In particular DW(h) — h for every h £ H. Notice that in this particular case we have an explicit rela- 
tion between the covariance of W and the inner product of the Malliavin derivate, Cov[W(h 1 )W(h 2 )] = 
EiW^W^)] = (hx, h 2 ) n = (DWih^DW^n. 

Remark 2. In the case of a centered stationary Gaussian process, X t , the Hilbert space can be chosen in 
the following way: 

Consider the inner product (l[o,t]> l[o,s])w = Cov[X t X s ] = C(t — s) and take the Hilbert space T-L as the 
closure of the set of step functions on K with respect to this inner product. With this Hilbert space one has 
that X t — W(l[o,t]) and DX t = ho,t] where D is the Malliavin derivative. 

Since the Malliavin derivative verifies the chain rule, we have Df(F) = f'(F)DF, for any / : R — ► R of 
class C 1 with bounded derivative (is also true for functions which are only a.e. differentiable, but with the 
assumption that F has absolutely continous law). The second Malliavin derivative (denoted by D 2 ) can be 
define recursively. For k > landp> 1, O k ' p denotes the closure of S with respect to the norm ||-|| fc defined 

by 

\\ni P = E[\ F \*>]+j2E[\\ Dl F\\U] 

i=l 

Consider now an orthonormal system in H denoted by {ej,/ k > 1}. Then, given elements <fi £ l-l® kl ,ip £ 
V. (S ' k2 , the contraction of order r < min{fci, k 2 } is the element of %®( fc i+ fc 2-2r) defined by 

oo 

i\ ,...,i r — l 

In particular, <fi <g) r %j) = ((f), ip)-H® r wnen k\ = k% = r. 

Remark 3. Again, in the white noise framework (when H = L 2 J, for symmetric functions 4> £ L 2 iS)kl , 
ip £ L 2 ®^, the contraction is given by the integration of the first r variables, i.e., (j)® r ijj = ((/>,ip) L 2 . Also, 

we have a formula for the product of stochastic integrals]^ I p (f)I q (g) — X)r=o r '(r) (r) Ip+q—2r (/ ®r <?) 

Define the divergence operator 8 as the adjoint of the operator D, so if F £DomS then S(F) £ L 2 (il) 
and E[8{F)G] = E[(DG, F) H }. 

Remark 4. When H = L 2 , the divergence operator 8 is called the Skorohod integral. It is a creation operator 
in the sense that for all F £Dom8 C L^ xP (T x f2) with chaos representation F(t) = X)^=o ^<?(/<?(^ ')) 
(/ 9 G i 2 i8(g+i) are symmetric functions in the last q variables), 8(F) = J^L ^<?+i(/<?)0 

For all F e L 2 (51) denote by J q F the projection of i 7, in the q th chaos. Then, the Ornstein-Uhlenbeck 
semigroup is the family of contraction operators {T t / t > 0} on L 2 (fi) defined by T t F = YlqLo e ~ 9t JqF- 
Using Mehler's formula we can find an equivalence between Mehler's semigroup and the O-U semigroup. 
More formally, take W' as an independent copy of W defining (W, W') on the product probability space 
(fl x VI, IF® F', P x P). Each F £ L 2 (Q) can be regarded as measurable map F(W) from R w to R determined 
Poff-ks. such that T t F = E'[F{e- t W + Vl - e~ 2t W')] . The infinitesimal generator for this semigroup 
(denoted by L) is given by LF = Y^Lo^lJqF and DomL = B 2 ' 2 =Dom(<5.D). It can be proved that, 
for F eDomi, 8DF = —LF. The pseudo-inverse of this operator (denoted by L^ 1 ) is given by L~ 1 F = 
T,T=i T" J 9 F ' and is such that L lF eDomi and LL ~ lp = F - E[F] for any F £ L 2 (n). 

1 a A b = min{a, b} and a V b = max{a, b} 

2 /is the symmetrizatoin of /, i.e., f{zi,. . . , z q ) = i J2„ f( z a(i),- ■ ■ ,Z*(q)) 
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The use of Hermite polynomials is extremely important in this setting. The relationship between Hcrmitc 
polynomials and Gaussian random variables is the following. Let Z\, Z 2 be two random variables with joint 
Gaussian distribution such that E[Z\\ = E[Z$] = 1 and E[Z X \ = E[Z 2 ] = 0. Then for all q,p > 0, 

E[H p (Z 1 )H g (Z 2 )] = \f E ^ q (3) 

On the other hand, it is possible to expand a C 2 function / : M — > K in terms of Hermite polynomials, that 
is, 

OC 

f(x) = E[f(Z)]+Y / c q H q (x) (4) 

-3=1 

where the real numbers c q are given by c q q\ = E[f(Z)H q (Z)] and Z ~ A/"(0, 1). 

Remark 5. Notice that in the white noise case we have the relationship H q (W(h)) — J Tq h^ q dWt 1 ■ ■ ■ dWt q — 
Iq(h® q ) so the decomposition of f(W(h)) can be regarded as 

oo 

f{W{h))=Y J C q I q { hm ) 

9=0 

With © and ((H) we are able to compute the covariance of a real function / in the following way, 

oo oo 

Cov[/(^)/(^ 2 )] = J2 ^ E [ H p{Zi)H q {Z 2 )] - ]T c 2 M E i Z i Z ^) 9 ( 5 ) 

p,q=l q=l 

With this background, is possible now to state the main tool developed in |20j and used in this paper. 
Lemma 1. Let F G ED 2 ' 4 with E[F] = [i and Var[F] = a 2 . Assume that N — W(/i, a 2 ), then 



dw(F, N) < ^E[\\D 2 F ®! D 2 F\\ 2 nm ] » x E]\\DF\\%\ * (6) 

where dw is the Wasserstein distance given by dw(U,N) — su P{ ff -|| ff || L . <i} \E[g(U)} — E[g(N)]\. If, in addi- 
tion, the law of F is absolutely continuous with respect to the Lebesgue measure, then 



drv(F, N) < y^-E[\\D 2 F ®! D 2 F\\ 2 Hm ] * x E[\\DF\\ 4 n ] * (7) 

where drv is the total variation distance given by drviU, N) = sup^ e g/ R N \P(U G A) — P{N G A)\ 

Remark 6. This lemma is quite useful since it essentially tells us that in order to obtain a CLT result for 
a family of random variables Ft , we just need to check three conditions: 

1. Expectation of the First Derivative's Norm: 

E^DFtW^} = 0(1) as T^oo (8) 

2. Expectation of the Contraction's Norm: 

E\\\D 2 F T ®i D 2 F T ||^ 02 ] -> asT^oo (9) 

3. Existence of the Variance 

Var[F T ] -> S 2 G (0, oo) exists as T — > oo (10) 



Due to the Gaussian Poincar'e inequality, Var[FT] < £'[||.D.Fr||?{] < y E^\\DFt\\^i\ > 80 ^ e var ^ ance w Hl 
go to if the expectation of the first Malliavin derivative 's norm goes to 0. This is why condition J5J) is 
necessary, and the convergence to zero of the Wasserstein distance relies on condition {Pp. 
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2.2 Malliavin calculus on Wiener-Poisson space 

Let T~L = Lfr. Assume there is a complete probability space P) where L t is a cadlag, centered, Levy 

process: L t has stationary and independent increments, is continuous in probability and Lq — 0, with 
E[L 2 } < oo. At the risk of causing some confusion, is denoted by T the filtration generated by L t completed 
with the null sets of the above filtration, and work on the space (f2, J 7 , P). Assume this process is represented 
by the triplet (0, c 2 ,^), where v is the Levy measure s.t. d/u,(t,x) = a 2 dtSo(x) + x 2 dtdv{x){l — 5o(x)) and 
f M x 2 dis(x) < oo. This process can be represented as 

L t = aW t + J J xdN{t, x) 

where Wt is a standard Brownian motion, a > and N is the compensated jump measure. A fuller exposition 
on Levy processes can be found in [2] and |33) . This process is extended to a random measure M, which is 
used to construct (in an analogous way to the Ito integral construction) an integral on the step functions, 
and then by linearity and continuity it is extended to L 2 (([0,T] x R) q , B([0, T] x K)«,/x®«) and denoted by 
I q . This integral satisfy the following properties: 

1. Iq(f) = IqCf) 

2. Iq(af + bg) = al q {f) + bl q (g) (a, b G E) 

3. E[I p (f)I q (g)] = q\J im+xM)q fgd^n {q=p} 

These properties are stated in jM] and their proof can be found in Q3]. We have a product formula in 
this framework which is similar to the one in remark [3] but with extra terms coming from the poisson 
integration part. A product formula for the pure jump framework can be found in [25]. Before stating 
the formula, is needed to define a general version of the contraction. Let <fi £ L 2 mkl and ip € i 2 1 «* !2 be 
symmetric functions. Then the general contraction of order r < min{fci,A;2} and s < min{fci,fc2} — r is 
given by the integration of the first r variables and the "sharing" of the following s variables, i.e., <fi ip = 
n*=i*2i ((f>(-,z,x),ip(;z,y)) L 2 , where z G (M 2 ) s and (x,y) G («. 2 ) kl - r - s x (R 2 ) fc2 ~ r - S . Now the product 

formula can be stated as follows. If |/| ^ \g\ G L 2 iP+g _ 2r _ s for < r < min{p, q} and < s < min{p, q} — r 7 
then 

pAqpAq-r / \ / \ / \ / \ 

w)im=j: e ri - i (;)(;)( p 7 )( 9 7 jw^/®^) w 

The proof of this product formula can be found in [17] . 

Remark 7. In the general contraction formula z-n is the size of the jump and zn-\ is the time when that 
jump occurs (z = (zi, Z2, ■ ■ ■ } Z2 S -i, ^2s)) ■ If we only have the Wiener part, the factor zn would be zero 
unless s — 0, and we then obtain the contraction defined in the Wiener space. Similarly, when the terms 
where s ^0 are zero, the formula ill)) reduces to that in remark\3[ 

It has been also proved by Ito [HJ Theorem 2] that for all F G L 2 (0) := i 2 (0, J 7 , P), we have, 

oo 

F = J2 I M> / ? e^,:=I 2 (([0,r]xM)',B([0,T]xM)«,^) (12) 

9=0 

and that this representation is unique if the / g 's are symmetric function. From this chaotic representation 
we can define the annihilation operators and creation operators, the former will be the Malliavin derivative 
and the latter will be the Skorohod integral. In this way define DomD as the set of functionals F G L 2 (il) 
represented as in (jT2j) such that Ylq*Li 11^- ll/Jlz^ < oo. For F GDomZ) the Malliavin derivative of F is the 
stochastic process given by 

oo 

D Z F = 2 g/ g -i (/,(*,-)), z G [0,T] x R, /, symmetric (13) 

? =0 
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If we define the inner product as (f,g) L 2 = j R+yM f(z)g(z)d^i(z), then DomD is a Hilbert space with the 
inner product (F,G) = E[FG] + E[(D Z F, D Z G) r2 }. We can embed DomD in two spaces DomD and 
DomD J . DomD is defined as the set of all functionals F G L 2 (Q) with representation given as in (|12 j) such 
that ^2qLi Jo 0)i ') H^a ^ < °°; while DomD' 7 is defined as the respective functionals satisfying 

Sq=i ??• /[otIxRo ll/«( z ' ^m( z ) < °°; hence DomD=DomD°nDomD J . We can now rewrite (due to 
the independency of 14 7 and iV) ft as the cross product flw x ^jv- 

• Derivative D^o 

This derivative can be interpreted as the derivative with respect to the Brownian motion part. Using 
the isometry L 2 (fl) ~ L 2 (Vlw] L 2 (il^)), we can define a Malliavin derivative as we did in the Wiener 
case but using the L 2 (f^)-valued smooth random variables Sfj, that is, for the functionals of the 
form F = Y^i=i GiHi, where Gi G S and Hi € L 2 (£l^). Then, this derivative will be D w F = 
'Y^i = i{D W Gi)Hi and this D w Gi is the derivative defined in ([2]). This definition is extended (see [34] ) 
to a subspace DomD 11 ' cDomD and for F EDomD w , 

Dt,oF = -DY F (14) 
a 

Furthermore, we also have a chain rule result for functionals of the form F = f(G,H) G L 2 (il) with 
G GDomD^ 7 , H G L 2 (fl^) and f(x, y) continously differentiable in the variable x with bounded partial 
derivative. We have that F GDomD and 

D tfi F=-^(G,H)DYG (15) 
a ox 

This is also true (as in the Wiener space case) for functions which are a.e. differentiable but with the 
restriction that G has an absolutely continuous law. 

• Derivative D z (z ^ (t,0)) 

This derivate has been shown, in |35| . to be a difference operator. The idea is to introduce a jump 
of size x at moment t. Then, the Malliavin derivative with z = (t, x) will be the translation operator 
given by 

x 



for any F eDomD" 7 such that -E[J/j T ^ xRo (^ z F) 2 d^i(z)] < 00. See [35] for a complete contruction on 
the canonical space in which this is developed, and [34] for a quick explanation on how to introduce 
a jump at moment t and the conditions on the cj's. We also have a chain rule for this Malliavin 
derivative but using just the difference instead of the derivative, that is, for F — f(G, H) g L 2 (Vl) with 
G G L 2 (Qyy), H GDomD" 7 and f(x,y) continous it holds that, 

/(G, H{bj t , x )) - f(G, H(u>)) f(G,xD t H + H(u))-f(G,H(u)) 
L) z b = = 

x x 

Notice that if / is differentiable, then we can use the mean value theorem to obtain 

D z F=^(G,H(uj) + e z xD z H)D z H (16) 

for some 9 Z G (0, 1) 

In the same way, consider the chaotic decomposition F(z) — Y^Lo Iq(fq( z i ')): with / G L 2 ^ %q symmetric 
with respect to the last n variables. If Y^Loil + -0' HAHl 2 < 00 then we say that F GDom<5. Now we can 
define the Skorohod integral of F GDom<5 by 

^)^j, +1 (/,)a 2 (si) (17) 

9=0 
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This operator is the adjoint of the operator D z , so E[5(F)G] = E[(F(z), D Z G) L2 ] for all G eDonxD. 
Denote by L 1 - 2 the set of elements F G L\ *®p([0, T] x M x ft) such that QQ* ll/gllz, 2 < 00 • For all 

F G L 1 ' 2 cDoni(5 we have that F(z) eDomD, Vz/i- a.e. and that -D.F(-) G £ 2 82xP (([0, T] x K) 2 x Q). 

Finally, the definitions of the Ornstein-Uhlenbeck semigroup T t and its infinitesimal generator L are the 
same as in the Wiener space case. Basically, all we need to define it is the Malliavin derivative and the 
Skorohod integral, that is, we can just define L = SD. With this definition we obtain that for F E L 2 (fl) 
with chaotic representation (|T2"1) . LF — J2 q % Li~1^q(fq) an d T t F = YlqLo e ~ qt Iq (/<?)• Samely, the pseudo- 
inverse is given by L~ X F = £~ x =±I q (f q ) and LL~ X F = F- E[F]. 



3 Main theorems 



The first tool needed is the extension of the so-called Gaussian Poincare inequality for the present context. 
But to prove this is required to have an inequality similar to the one proved in [201 Proposition 3.1] (was 
proved for all p > 2 in the Wiener space case) . The technique used in their proof was based on the equivalence 
between Mehler and Ornstein-Uhlenbeck semigroups for the Gaussian case, but in the Wiener-Poisson space 
we lack such an equivalence. Nevertheless, it is possible to prove it for p — 2 and that is, in fact, the one 
needed to prove the extension of the Gaussian Poincare inequality. 

Proposition 1. Let F EDomD such that E[F] = 0. Then, 

E[^DL- l F\\ 2 L2 ] <E[\\DFf L2 ] 
Proof. Assume F has its chaos decomposition given by (fT2|) . By the orthogonality between chaoses we get, 

2 



sOlDz-^iy =E 



9=1 q L 2 ] l q=i q 



'■?;Jli= 



< E 



Y,\\ D Wq)\\ 2 Ll 

9=1 



E[\\DF 



□ 



(18) 



Theorem 1. (Extension of the Gaussian Poincare inequality) 
Let F EDomD. Then, 

Var[F] < E[\\DF\\ 2 L2 ] 

with equality if and only if F is a linear combination of elements in the first and th chaos. 
Proof. Assume, without loss of generality, that E[F] =0. 

Var[F] =E[F 2 ] = E [(DF, DL~ X F) L% ] < E[\\DF\\ 2 L2 ] ^E[\\DL,- l F\\ \ 2 ] ' < E[\\DF\\ 2 L2 ] 

where Proposition [T] was used in the last step, and the fact that F = SDL~ 1 F in the second step. □ 

Notice is possible to combine formulas (|15|) and a version of (|16|) to write the chain rule in a unique way 
when / is a function of just one variable which is k times continously differentiable. Similarly, is possible to 
give a unified formula for the derivative of a product. Using the fact that the jump x is zero for Dt,o, the 
following unified chain and product rules are obtained, 

Proposition 2. (Chain rules and Product rule) 

Let F, G, H £DomD w C\DomD J such that DF,DG,DH G L 2 . Also consider f,g G C fe ~\ both with bounded 



8 



first derivatives (f and g' can be unbounded if the the law of F is absolutely continuous with respect to the 
Lebesgue measure), such that f is k-times differ entiable and g( k ~ l > is a.e. differentiate. Then, 



(19) 



n—1 

for some function 6 Z G (0, 1) for all z = (t, x) £ R + x M, 

DMF) _ g £^- 1(D , Fr + £l!^l^.-. (Dlf _ (20) 

„ =1 ™ Jo l fe 1 J- 

and 

D Z (GH) = D Z G H + G- D Z H + x ■ D Z G ■ D Z H (21) 

for all z = (t, x) G M.+ X M. 

Proof. If z = (t, 0), we get the chain rule formula (fT5|) and the usual product rule in Wiener space. If 
z — (t, x) with x 7^ 0, then the Malliavin derivative formula tells us that D z f(F) = ^ F ^ t - :c ^~^ F ^ u ^ , Since 

/ g C fe - X and fc-times differentiable, f(y) = f(y ) + ^^(y - 2/o) n + - y ) k for 

some € (0, 1) (Taylor series of / with mean-value form for the remainder). Using this expansion with 
y = F(u> z ), yo — F(uj) and recalling that y — yo — F(lu z ) — F(lj) — xD z F, the chain rule result (fT9|) will 
follows. 

The second chain rule formula is obtained by using the Taylor expansion for g with integral form for the 

remainder, i.e. g(y) = g(y Q ) + Y!ZX 9 J (V ~ + fy ffc^Tn ~ v) fc_1 cfc- By using the values for y 
and yo as in the previous case, and applying the change of variable v = F + xu in the integral we get the 
chain rule formula (|2"U|) . 

For the product rule we get trivially that 

DJGH) = G (^,,)H(jo t , x )-G(u)H(uj) = G(u t , m )-G(.u) H ^ + )-H(u) + x G(u t . x )-G(u) Hj^-Hju) 

□ 

3.1 Upper bound and second order Poincare inequality 

The main theorem of this paper is now stated. Recall the following bound on the Wasserstein distanc^l, 

d w (F,N) < sup \E[f'(F) — Ff(F)]\ 

where N ~ Af(0, 1) and J w := {/ G C 1 / /' is Lipschitz, Wf]^ < 1, \\f" ^ < 20 
Theorem 2. (Upper Bound) 

Let N ~ 7V(0, 1) and let F eDomD w C\DomD J be such that E[F] = 0. Then, 

d w (F,N)<E^l-(DF,-DL- 1 F) L2 + E[(\x(DF) 2 \ ^DL' 1 F\) L2 ] (22) 

Proof. By Proposition [2] (with k = 2 in flU) we get Df(F) = f'(F)DF + f" F f"(F + xu)x(DF - u)du. 
On the other hand, using the identity F = — 8DL^ X F (recall E[F] = 0) and the integration by parts formula 



3 see 1191 for further details on this bound 

4 The class Jyy is the collection of all continously differentiable functions / : R — > R that has derivative bounded by 1 and 



such that there exists a version of /" that is bounded by 2 
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we get E[Ff(F)) = E[(Df(F), -DL^F) J. Putting them together we get 
\E[f'{F) - Ff(F)]\ = E[f'(F)] - E[(Df(F), -DL^F)^] 



E[f(F)} - E 



/ f DF \ 

/'(F) (DF, -DL^F)^ + (J f"(F + xu)x{DF - u)du, -DL^F 



< E 



\f'(F)\ l-(DF,-DL- 1 F) 



Li 



<Wf'\\ooE 



1 - (DF,-DL- 1 F) 1 



E 



E 



/ \f"{F + xu)\ \x(DF — u)\ du, \DL~ 1 F\ 

„,„,, HD z Ff\ 



\DL~ F\ 



L'i 



Finally, use the fact that H/'H^ < 1 and ||/"|loo < 2 to obtain the result. 



□ 



Is desirable to use this result to obtain a nice upper bound, as in the Wiener space with Lemma [TJ The 
main problem faced is that equivalence between the O-U and Mehler semigroups is no longer available, so the 
proofs of FlIlD 2 !,- 1 ^!^] < ±E[\\D 2 F\\ P L2 ] and E[\\D 2 L' 1 F\\ P op ] < ±E[\\D 2 F\\ P Qp \ given in the Wiener 

space case fail for the Wiener-Poisson space case. Nevertheless, is still possible to state an equivalent version 
of Lemma [T] for the case when F lies in one specific chaos. 

Corollary 1. (Second order Poincar'e inequality) 

Fix <? G N and let F = I q (f) with E[F] = \i and Var[F] = a 2 . Assume that N — Af(fJ,, <J 2 ), then 



d w (F,N)<^(2E^I%I^E[ DF\\l^ E 



(x,(D 2 F) 2 ) 



Li 



-^E[(\ X \,\DFf 



LV 



(23) 



< ^[2E[\\D 2 F(g> 1 D 2 F\\' L2 _] i E[\\DF\\ i Ll ] i +E 



(x,{D 2 F) 2 ) 



Li 



qa° 



E\( \x\ , \DF\ 

(24) 



Proof. Assume without loss of generality that /i = and a 2 = 1 . By Theorem [2] and Holder we have that 



dw(F,N) <E[{1- (DF, -DL^F)^] 5 + E[(\x(DF) 2 \ , {DL^F])^] 

Also, notice that E[(DF,-DL- 1 F) r2 ] = E\—5DL~ 1 F ■ F] = E\F 2 } = 1, so if G = (DF,-DL- l F) r2 

then E[(l - (DF, -DL^F) 2 ) 2 ] =Var[G]. By Theorem[I]we have that Var[G] < E[\\DG\\ 2 L 2 ]. Also, by 
the product rule (f2"Tl) we have that 



DG = (D 2 F, -DL- X F) L2 + (DF, -D 2 L- 1 F) L2 + (xD 2 F, -D 2 L- 1 F) L2 
Putting all together and using the fact that —L~ 1 F = ~F we get, 

E[(l - (DF, -DL~ 1 F) L2 ) 2 ] ^ <—(2e\ (D 2 F,DF) 



1 \ 



Li 



E 



(x,(D 2 F) 2 ) 



Li 



The first term in the right is bounded above in the following way 



E 



(D 2 F, DF) 



Li 



<E[\\D 2 F\\J*E[\\DF\\%V 
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since 



(D 2 F,DF) L2 < \\D 2 F\\ o \\DF\\ 2 L 2 and by Holder. On the other hand, and again using the fact 
that F is in the q th chaos, we get 

E[(\x(DF) 2 \, \DL- l F\) L2 ] = -E[(\x\ , \DF\ 3 ) ] 

thus, obtaining the first inequality. For the second inequality, it suffices to see that if {jj}j>i is the sequence 
of random eigenvalues of the random Hilbert-Schmidt operator / — > (/, D 2 F) , 2 , then 



|D 2 F||^ p = max | 7 / < £ | 7 / = \\D 2 F 01 D 2 F 

3>1 



□ 



As in Lemma [T] this corollary basically says that if we want to show a CLT for a family of random 
variables Ft (living in a fixed chaos) it is sufficient to check the following conditions, 



1. Expectation of the First Derivative's Norm: 

B[||£)F r ||* 2 ] =0(1) as T->oo 

2. Expectation of the Cube of the First Derivative's Norm: 

E[(\x\,\DF T \ 3 ) ]-*-0 as T^oo 



(25) 



(26) 



3. Expectation of the Contraction's Norm 

E 



D 2 F T (g> 1 D 2 F T \\ 2 L2 



— ¥ as T^oo 



4. Expectation of the Squared Second Derivative's Norm: 

2 



E 



(x,{D 2 F T ) 2 ) 



LI 



L'i 



— > as T -> oo 



5. Existence of the Variance 



Var[F T ] — > £ 2 <E (0, oo) exists as T ->• oo 



(27) 



(28) 



(29) 



4 Special cases and applications 

4.1 The Wiener space case: 

Linear functionals of Gaussian-subordinated fields 

When we are working in this space the jump size is always zero, so the upper bound for the Wasserstein 
distance becomes 



d w (F,N)<E I - (DF.-DL^F) 2 
which coincides perfectly with the bound computed in [19] . 



(30) 
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Remark 8. It is important to stress that Theorem\^is not a direct extension of the inequality in [19], even 
though they appear to be similar. This is because in the pure Wiener case, the Malliavin calculus theory is 
developed for more abstract Hilbert spaces than L 2 , and the inequality proved therein, 

d w (F,N) < E\l- {DF,-DL- l F) n \ 

holds for any Hilbert space H. However, when H — LZ, Theorem^ is indeed an extension of Theorem 3.1 
in EE/. 

The Wiener space is well understood and the upper bounds obtained from (|30[) are more powerful than 
Corollary [TJ In fact, Lemma [TJ is true for all F £ O 2,4 and not just for functionals in a fixed Wiener chaos. 
As an application of Lemma [TJ the authors of |20) proved a very useful central limit theorem for linear 
functionals of Gaussian-subordinated fields. Before stating it, some notation is introduced: Let X t be a 
centered Gaussian stationary process and define C(t) = E[XoX t ] = E[X s X t + s ], its covariance function. By 
remark [21 is known that the Malliavin derivative of X t is well defined. Let T > 0, Z ~ A/"(0, C(0)) and 
/ : R -> R be a real function of class C 2 not constant s.t. < oo and i?[|/"(Z)| 4 ] < oo. In order 

to simplify the notation, the following random sequence is defined, 

F T =T"5 [ T (f(X t )-E[f(Z)])dt 
Jo 

The theorem is stated as follows, 

Lemma 2. Suppose that L |C(i)| < oo, and assume that f is asymmetric real function. ThenliuiT-^oo Var[Fx] := 
E 2 £ (0, oo) exists and as T — > oo 

F T ^N~M(0,Z 2 ) 

Remark 9. One advantage of inequality {6p is the fact that it allows us to quantify rates of convergence to 
normality. Indeed, it has been proved (see \20/ ) that if Z ~ A/"(0, 1), then (as T — > oo), 

* r (7^' , )- ,r,) 

Our goal in this subsection is to extend this result to the case when L |C(t)| = oo. This is achievable un- 
der some conditions on the decay rate of the covariance. In fact, it is very convenient that for this functional 
the conditions JSJ), © and (fTUt reduce to just one condition on the covariance of the underlying stationary 
Gaussian process X t . Let V(T) be a strictly positive continuous function with V(T) — > as T — > oo such 
that cither TV(T) -> or V £ C 1 and TV'(T) -> as T -> oo. The following is the condition on the 
covariance that replace the three conditions on remark [51 



Condition *: Either J R \C(t)\ < oo or V(T) (with the above characteristics) exists such that, 

V(T) T^oo r 

V(T) represents the decay rate for the covariance function. Consider the following function 



V(T) 



'T if f™\C(x)\dx < oo 

Jo lo V(x)dxdy if f™\C(x)\dx = oo 



Let M c := {/ £ C 2 / f is symmetric if J R \C(t)\ < oo or E[f(Z)Z] ^ if J R \C{t)\ = oo} and rewrite the 
functional Ft as follows, 

F T = V(T)-i [ (f(X t )-E[f(Z)])dt 
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Theorem 3. Suppose that condition * is verified by C(t) and that f G Mc- Then lim-r-s-oo Var[.Fr] := S 2 e 
(0, oo) exists and as T — > oo 

F T ^> N ~ Af(0,Z 2 ) 
Furthermore, if J R \C(t)\ = oo, then S 2 = 2M (E[f{Z)Z]) 2 . 

Before tackling this theorem, it is necessary to verify some facts that would simplify the proof. 
Proposition 3. Suppose that J R \C(t)\ = oo. Then as T — > oo, 

1. (fiv(x)dx)- 1 fi\C(t)\dt = 0(l) 

2. y(T)- 1 J [0T]2 \C(t - s)\dsdt = 0(1) 

3. • IfTV(T) -> 0: 

V(T)- 2 T(fi \C(t)\dtf = 0(m ax {V(T),TV(Tr(^V(x)dx)- 1 }) 
• IfTV(T) -» and TV'(T) -> 

y(T)- 2 T(/ T |C(t)|dt) 3 = 0(max{V(T),7V(T)}) 

4- For fixed q > 1: 

F(T)- 1 / C(i - *)«d«tt 2M1 {?=1} = T M # 9 = | 
J[o,iY (0 iff q^l 

Proof. The proof just involves simple applications of L'Hopital's rule (L). 

Bm £m>\* , lim jcmi _ |M| 

T ^°° Jo T— >oo V(T) 

2. Notice first that Jj T j 2 |C(t — s)| dsdt = 2 f Q T |C(a;)| dxdt so 

um W ic (( -.)i^ _ um 2^j'|cmi^» j. 2 Bm jcp _ 

r-^oo y(T) ^oc J q t fV V {x)dxdy T ^ V(T) 

3. • If 7V(T) 0: 



=0(1) 

• If TV(T) and TV'(T) -> 



T(fn\C{t)\dt) 3 ( \l '\C{t)\dt\ 3 T(^V{x)dxf l , , N VN 

lim Uo J Wl ; = lim J ° T ' V n lim Uo „ V ; ; = 0(1) lim (W(T)+3TV'(T)) 



V{TY T^°°\ V(x)dx J T ^°° V(T) 2 T— >oo 



=0(1) 



4. If for q > 1, cither liniT^oo Jj T p C(i — s) q dsdt < oo or lim^^co Jq T C{x) q dx < oo, then the result 
will follow trivially. So let's assume that both go to infinity as T goes to infinity. 

L T]2 C(t- s)«dsdt 2 f n T fiC(x)«dxdt L C7(T)« fC(T)\ q , w n 

Um J|0 ' T| „ = lim J ° Jo V ; = 2 lim = 2 lim -i-i F r ^ 1 ) = 2M1 { ,_ 1) 

y(T) t^cx> fTfy V (x)dxdy r-oo V(T t^oo\V(T)J^ . ^ 



->0 if q>l 
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□ 



Now, to the proof of Theorem [3] 

Proof. Notice that if L \C(t)\ < oo then Theorem [3] reduces to Lemma [5] and there is nothing left to prove. 
Assume then, that L |C(f)| = oo. Due to remark [5J it is enough to check that condition * implies conditions 
©, © and (HUD- 

• Expectation of the First Derivative's Norm: 

— First Malliavin Derivative: 




— Norm of the First Malliavin Derivative: 



\\DF T \\ 2 n = ViT)- 1 [f(X t )f'(X s )(l m ,l [0 , s] ) n dtds = V^)' 1 [f(X t )f(X s )C(t - s)dtds 




Then. 



\\DF T f n = V(T)- 2 [f(X t )f'(X s )f'(X u )f'(X v )C(t-s)C(u-v)dtdsdudv 



J[0,T] 4 



Expectation of the First Malliavin Derivative's Norm: 

By using Holder (twice) on the expectation and by the stationarity of X t we have the bound, 



\E[f(X t )f'(X s )f'(X u )f'(X v )]\ < E[\f'(Z)t] 



finally recovering the power we get, 




— O(l) by Proposition [3] 



All this proves that, 



£[||DF r ||^] * = 0(1) as T->oo 



• Expectation of the Contraction's Norm: 

In the same way we get, 



— Second Malliavin Derivative: 




— Contraction of Order 1: 



D 2 F T ®! D 2 Ft =V{T)- x f f"(X t )f"(X s )l m ® l [0i>] (l [0|t] , lp,,]) dtds 



J[0,T] 2 



^V(T)- 1 [f"(X t )f(X s )l [0tt] ®l [0tS] C(t-s)dtds 



J[0,T] 2 
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Norm of the Contraction: 



\D Ft ®i D Ft\ 



V{T)- 2 f"(X t )f"(X s )f"(X u )f"(X v )C(t - s)C(u - v) 

J[0,T] 4 

x (l[o,t], 1 [o,u])n (Mo,a]A[o,v])n dtdsd ' udv 
V(T)- 2 [f"(X t )f"(X s )f"(X u )f"(X v )C(t - s)C(u - v)C(t - u)C{s - v)dtdsdudv 

J[0,T] 4 



Expectation of the Contraction's Norm: 

By using Holder in the same way as above we get, 

E[\\D 2 F T ®i D 2 F T \\ 2 2 ] < £[|/"(Z)| 4 ]i/(T)- 2 / \C(t - s)C(u - v)C(t - u)C{s - v)\ dtdsdudv 

./[O.T] 4 

Now, let's make the changes of variable y = (t — s,u — v,t — u,v), and let's denote the new region 
by h x [0,T]. So, 

BfllD 2 ^®!^ 2 ^!!" ] <i?[|/"(^)| 4 ]nT)- 2 / T / |C( 2 /i)C(y 2 )C(y3)C(y 2 + y3-yi)|rf2/ 

= E[\f'(Z)\ 4 ]V(T)- 2 T [ \C(y 1 )C(y 2 )C(y 3 )C(y 2 + y 3 - m )\ d Vl dy 2 dy 3 
Jn 

Taking in account that by Cauchy-Schwarz, V t e K, 

<i 



=C(0) 



\C(t)\ = Jg^L ^Var^olVar^] < C(0) 
VVar[X ]Var[X t ] 

Also, it is clear that C [— T, T] 3 and since the integrand is a non-negative even function, we can 
deduce that, 



E\\\D 2 F T ® l D 2 F T \\> um ] <E[\f{Z)\*]C{V)V{T)- 2 T[2 ] \C(v)\dy 

= 8E[\f"(Z)\ 4 ]C(0) V(T)- 2 t( f T \C(y)\ dy ' 



-^0 by Proposition [3l 



All this proves that, 
* If TV(T) ->■ 0: 



£[||£> 2 F T 1 .D 2 J F r ||J i82 



* = O ^max I V(T), IV(T) 2 ^ V(x)daA 



as T — > oo 



If TV(T) and TV'(T) -> 



EfH^T 0i -D'^t ||^ 82 ] 4 =0(max{y(T),TF'(T)}i) as T 



• Existence of the Variance: 

Since / e A4c then E[f(X )X ] = E[f(Z)Z] ^ 0. Also Hi(x) = x, so the first Hermite constant in 
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the expansion (|4} is not 0, i.e., ci = E[f(Xo)Xo] ^ 0. Using the formula (0 for the covariance of / 
we get, 



Var[F T ] = E 



V(T)~iJ o {f(X t )-E[f(Z)})dtj 
ViT)- 1 ^c 2 q q\{E[X t X s \) q dtds = ^clqW{T)- 1 C(t-s) q dtds 

J[o,t] 2 q=1 q=1 J[o,r\* 

a L T p C(* - S)dsdt ~ f T] 2 C(* " 

Cn — — - — h > c„o! — — ■ — > 2Afc 1 

V{T) ^ 9 U(T) 



V(T)" 1 / Cov[/(X t )/(X s )]dids 

J[0,T] 2 



— ¥2M by Proposition [3] — >0, V q by Proposition l3l 

All this proves that, 



lim Var[F T ] = 2M (E[f(Z)Z]f e (0,oo) exists 



Since the conditions were satisfied, Theorem [3] is proved. □ 

Remark 10. Notice that during the proof of this theorem was possible to establish an estimate for the 
convergence rate to normality like in remark^ i.e., if Z ~ -*V(0, 1) then as T — > 00, 

• If TV(T) -» 0: 

1 

oT \ -1] 4 



Ft 



VVarpV 
• If TV{T) and TV'(T) 0: 

4.1.1 Examples 



Z^j =0 ^max i V(T), TV(T) 2 (J V{x)dx 



d w _^_,Z =0(max{F(T),TF'(T)}i) 



According to Theorem[3]the only condition we need to check in order to apply the central limit theorem to Ft 
is the decay rate of the covariance function for the underlying stationary Gaussian process X t (condition *). 
In fact, if the decay rate is t~ a then we can apply the CLT if a € (0, 1) U (1, 2), because in the case a £ (1, 2) 
the integral L C(t) is finite and in the case a £ (0, 1) the same integral is infinite but V(T) = T~ a £ C 1 
and TV'(T) = -aT- a — > as T — > 00. 

1. Fractional Brownian Motion (FBM): 

In this particular case is known that the difference of FBM is a stationary process for all H £ (0, 1). 
So, if Bf is FBM then X t = B^ +l — Bf is a centered Gaussian stationary process. Its covariance 
function is, 

C X {T) = E[X T X ] = E[(B» +1 - B?){B? - <)] = 1 

Thus, 

(1 + JA 2H + (1 _ L) 2H - 9 
lim T 2 - 2H C\(T) = lim j' 2 A___ZIz + ^ ^ £ = ff(2H _ 1) = M G (0,oo) 

Then, the decay rate of its covariance function is t 2H ~ 2 , so Theorem [2] is applicable to the increments 
of FBM for all H £ (0, ±) U (§, l), and F T iV ~ 7V(0, E 2 ) as T -> 00 



16 



2. Ornstein-Uhlenbeck Driven by FBM: 

This process is given as the solution of the following SDE: 



Y t H =Y H -\fY s H ds + oB? 
Jo 



(31) 



where a, A > are constants, and Bp is a fractional Brownian motion with Hurst parameter H £ 
(0, |) U (4, l). This process is stationary due to the stationarity of the increments of the FBM (used 
in the first example). So X t — Y t H — E\Y^] is a centered Gaussian stationary process. In [5], the 
authors proved the following lemma, 

Lemma 3. Let H e (0, ~) U (|, l) and N € N. TTien as T -> oo. 



~ 2 /2n— 1 v 



n=l 



fc=0 



which basically tells us that for all H 6 (0, |) U (i, l) the decay rate of C^T) is very similar to the 
decay rate of C\(T) (the covariance of the FBM increments). Lemma [3] implies that, 

lim T 2 - 2ff C 2 (T) = H{2H ~ 1)?2 =M € (O.oo) 

T-yoo A z 

As in example 1, due to this rate of decreasing, Theorem [2] is applicable to this process for all Ft £ 
(0,|) U and F T ^ JV - A/"(0, £ 2 ) as T -> oo 

According to remarks |H] and [TU] we can tell that for the above examples Ft has a rate of convergence to 
normality of at least for all ^ £ ( 0) 1) {J ^ ^ that ig) for ^ _ ^ 



Ft 



VVar[F T ] ' 



,z) = o(t<- 



as T — > oo 



4.2 The Poisson space case: 

Simulation of small jumps 

Notice that in this space the measure [i has support contained in R + x M , and the formula obtained from 
the general case will be 



d w (F,N)<E \-(DF,-DL- 1 F) L2 + E[{\x(DF) 2 \ ^DL^Fl) L2 ] 



(32) 



Remark 11. This particular case was worked out by G. Peccati, J. L. Sole, M. S. Taqqu and F. Utzet in 
J29}j , with several examples (and conditions) given. They obtained the following inequality, 



d W {F,N) < E 



1 



DF, —DL~ 1 F 



E[{ (DF)' 



DL~ F 



(33) 



where the factor x of the second term is missing. It is important to stress that both formulas are equivalent 
and that the difference lies in the definition of the Malliavin derivative and the random measure. In fact, 
the definition of Malliavin derivative used by them is D Z F = F(w 2 ) — F(u>) and the random measure is 
J A dN(t,x) (instead of J A xdN(t,x) as in our case). So, the kernels of both chaos decomposition will differ 
by a factor of x, that is, if F has chaos representation as in \12) for our framework with kernels f q , then 



F = Y J I q {*f q ) 

9=0 
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in their framework (I q is the corresponding extension with their random measure). For example, the Ornstein- 
Uhlenbeck process (with inital value 0) is given by 



- x ^xdN(s,x) = / 1 ( e - A (*- s )) =I 1 (xe- x( - t -^) 



and it is easy to check that DY t = xe A ^ s ) = xDY t . Also, since dfj,(t,x) = x 2 d/i(t,x), we have that 



DF,-DL~ l F) , = {DF,-DL- l F) 



Li 



Li 



and 



{DFf 



DL~ Y F 



= (\x(DF) 2 \ , \DL~ 1 F\) 



l 2 , 



Remark 12. In ISO! the authors accomplish the remarkable generalization of this bound to the multi- 
dimensional case. Their theorem reads as follows, 

Lemma 4. Fix d > 2 and let C = {C(i, j)}a<i,j<d be a d x d positive definite matrix. Suppose that 
N ~ jV(j(0, C) and that F — [Fi, . . . , Fd) is a Mr -valued random vector such that E[Fi] — and F{ GDomD 
for all i = 1, . . . , d. TTienJf) 



d 2 (F,N)<\\c- 1 \\ op \\c\\iE 



Y,(C(i,j)-{DF i) -DL- 1 F j 



2tt 



^2\df ^Idl^f 



From this inequality they conclude a useful result for the first chaos case. 
Lemma 5. For a fixed d > 2, let N ~ A/d(0, C), with C positive definite, and let 

F„ = (F n ,i, . . .,F n ,d) = (Ji(/i n ,i),. . . ,h(h n; d)), n>l 
Let K n be the covariance matrix of F n , that is, K n (i,j) = E\l\ {h n ,i) I\ = \h n ,i,h ni j\ . Then 



Li 



-> for all i,j = 1 . . . d. 



F n ^NifK n (i,j) >C(i ,j) ant 



By remark we know that it is possible to rewrite these conditions in terms of "our" framework. 



Indeed, denoting h n ^ — xh n ^, it follows that K n (i, j) = (h nj i,h n j) = (h nj i,h n j) L2 and 



M 5 IVzl 5 



In PQ, the authors proved that the small jumps from a Levy process can be approximated by Brownian 
motion. Before this theorem is stated, some notation needs to be introduced: Let X t be a Levy process 
with triplet (b,o~ 2 ,v). To isolate the small jumps, consider the variance c(e) 2 = Jn x i< e \ x dv(x) and the 

small jumps process XI = o{e)~ x ff, Q t ] x {| x |< £ } xdN(s, x). So X t = b e t + oW t + N£ + a(e)X^ where N£ — 
t^2s<t ^^"sl{|AX,|>£> is the part of (finitely many) jumps bigger than e. Their theorem reads as follows, 

Lemma 6. X\ ^—t Wt (Brownian motion independent ofWt and N£ for all e) as e — > if and only if for 
each k > 0, a(K,cr(e) Ae) ~ c(e) as e — > 



3 see 30 for the definitions 
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They also proved that the condition er(fter(e) Ae) ~ <r(e) as e — > is implied by lim e £iM = oo. The 

importance of this lemma is that X t ~ b e t + a 2 + cr(e) 2 Wt + iV t c (for e small enough), and the latter is 
quite easy to simulate. 

The objective of this subsection is to extend this kind of result to functionals that are not necessar- 
ily Levy. To focus just on the jump part, let's assume, without loss of generality, that the triplet of 
the Levy process X t is (0,0, v). Then, it can be written as X t — ii(l[o,t]) = jji t ] xR xdN(s, x). De- 



fine a t (e) 2 = ||M[o,t]x{|x|<e}| 



jfjT-Q t ] x{ia:|<e} h t (s, x) 2 x 2 dv(x)ds for some h t £ L 2 ^ and consider the 



processes X t = 2i(/j t l[o,t]), = h(&t(e) 1 h t l[ ,t]x{\x\<e}) ■ Similarly, define a(e) 2 = / { | i|<t( i 2 rfi/(i) 
and XI = /i(o : (e)- 1 /i t l [0it]x{ | :! ,| <e} ). This means that X t = N f t + a(e)X% where = h{h t l[o,t]x{\x\>e}) 
has finitely many jumps. Let W be an isonormal Gaussian process with covariance structure given by 
E[W(f)W(g)] = f w+ f(s)g(s)ds. 



Theorem 4. For a fixed t, X\ Z ~ W(0, 1) as e -> if 

If[o,t] x {\x\<e}\ xh t( s > x )\ 3d »( x ) ds 



? t (e)3 

law 



e->0 



-> 



(34) 



Moreover, suppose that h t (s,x) = h t (s). Then X\ W(h t ) if 

I{\x\<e}\ X \ 3 



e->0 



> 



Proof. To prove the first statement, it is enough to verify that the upper bound 
to zero. Notice that 

^ 2 II _ i || 2 

II-D^IIl 2 = || CT «( e )~ ' l * 1 [0,t]x{|a:|<e}|| i 2 =1 



(35) 

goes to zero as e goes 



so E 



1 - (DXj, -DL-^Xt) Ll 



1 - WDX. 



t \\Ll 



E 



DL~ XI 



= E 



x\, 



= 0. On the other hand, 

iT[0,t]x{|iE|<e} 



DXl 



xht{s,x)\ du(x)ds 



a t (ef 



So bound (|32|) goes to zero as e goes to zero if (|34|) is true. To prove the second statement, it is necessary to 

prove that for any times {t±, . . . ,td} the random vector (X^ , . . . , Xf d ) —> (W(h tl ), . . . ,W(h td )). This is 
very easy to achieve thanks to Lemma [3] from remark [T^l because only two conditions need to be checked. 
Notice that for all e 

K e (i,j) =E[XlXt] = (o{ey 1 h u l[ 0M }y.{\x\<e}iG{ty 1 h t] l[ ,t ] ]y.{\x\<e}) T 2 

f 

= (/H 1 l[o,t 4 ],^ i l[o,^]> £a =^[W('H 1 l[o,t 1 ])W ? (^l[o,t,])] = C(i,j) 
so trivially K e (i,j) — > C(i,j) as e — ¥ oo. And for the second condition we have that 

I{\x\<e} \xfdv{x) 



Id* \a(e)- x htAu 



2 

L'i 



\hs\ ds- 



l ti 1 [0,t t ]x{\x\<e}\ 

"J, JO 

->0 by J35} 

Since both conditions were fulfilled, the conclusion of the theorem follows. 







□ 
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Remark 13. Notice that if ht(s,x) = ht(s) then 

vt(e) 2 =[[ h t {s) 2 x 2 dv{x)ds= ( h t (s) 2 ds f x 2 dv{x) = f h t (s) 2 dsa(e) 2 

J J[0,t]x{\x\<e} JO J{\x\<e} JO 

so it immediately follows that \3J$ is true if and only if \35)) is true. Assuming then, that i35\) is true, it may 

appear by the first part of Theorem^ that we can conclude the second statement, i.e., ^J*ht(s) 2 ds) 2 XI — 

XI — >■ Z, ~ 7V(0, 1) as e -t which is equivalent to X\ — >■ Z — 7V(0, J* h t (s) 2 ds) . Also, since E[X^Xl] = 

£ , [W / (/itl[o j t])M / (/i s l[o. s ])] for all e, X[ and W(htl[o,t]) coincide in their covariance structures. But this is 

not enough to prove that X\ W(ht); it only proves that, marginally (and for fixed t), X\ converges to a 
normal random variable, which does not imply that the joint distribution for all times is normally distributed. 
Then the conclusion does not follows from there. Hence the neccessity of Lemma\5[ 

4.2.1 Example: Fractional Levy Process 

Condition (1351) is quite easy to be verified. In fact, the measure dv(x) — \x\ ^ 2+ ^ lr a <a;<b\dx for <5 G (—1, 1) 
and a, b > (no jumps bigger than b or smaller than —a) is such that (|3"5|) is fulfilled. To check this, note 

that f {]xl<e} \x\ 3 dis(x) = ^ and a(e) 3 = (ff^) 1 , so a(e)- 3 J {]x{<e} \x\ 3 dv{x) = O(e^) and © is 
true. So, for the example, assume that the measure v is such that (|35p holds. 

1. Fractional Levy Process (FLP): 

There are two ways to represent a fractional Brownian motion as an integral of a kernel with respect to 
Brownian motion, and both deliver the same process (see [15] for a thorough explanation). One is the 
so-called Mandelbrot- Van Ness representation which is an integral over the whole real line with respect 
to a two sided Brownian motion. The other, is the so-called Molchan-Golosov representation which is an 
integral over a compact interval. In the Levy case, it is proved in |37j that these representations delivers 
different processes with very different characteristics. Because of this "non-uniqueness", the FLP 
generated by the Mandelbrot- Van Ness representation is called FLPMvN, and the one generated by 
the Molchan-Golosov representation is called FLPMG. It is known that FLP's have the same covariance 
structure as FBM. The advantage of FLPMvN over FLPMG is that the former is stationary and the 
latter is not (in general), as is shown in [37]. Nevertheless, since FLPMG is derived on a compact 
interval, Malliavin calculus can be applied to it. 

Consider X t as an FLPMG, that is, X t = h(Kf) where 

(K?,K?) L2 = \(\t\ 2H + \s\ 2H -\t-s\ 2H ) 

According to Theorem g] since §Q is true and h t (s,x) = K?(s), it follows that X\ W(K^) 
as e — > oo. But W(KH) — B^ is a fractional Brownian motion. In this case, we conclude that in 
order to simulate the paths of an FLPMG X t , we just need to fix e small enough, and simulate the 
finitely many jumps part N£ — Ix{K^ln x \ >e \) along with an (independent) FBM part Bf — W{Kf), 

laW _ , % rj 

because X t ~ + o{e)Bf . 

4.3 The Wiener-Poisson space case: 
Product of O-U processes 

Finally, the second order Poincare inequality developed in the combined space is used to obtain a CLT for 
mixed processes. First, notice that if we have a double Ornstein-Uhlenbeck (O-U) process as a sum of a 
Wiener O-U process Y t plus a Poisson O-U process Z t (independent of Y t ), it can be proved (in two different 
ways) that the functional Ft = J Q Y t + Z t dt converges to a normal random variable as T — > oo. The 
first way is to seperate Ft to the terms f„ T Y t dt and T - ^ jl" Z t dt and use the inequalities (|3"0"]) and (|3"!?]) 
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respectively to prove that each part goes to a normal. The other method is to use inequality (|22|) and just 
do one computation. The second way is clearly faster (since the kernels are the same). This is one advantage 
of having the inequality in the combined space. 

For our example, lets focus on a process for which we cannot use either of the inequalities, (|30l) nor (pl2"j) . 
to prove a CLT. First, for simplicity's sake (to avoid dealing with constants), assume that the triplet for the 
underlying Levy process is given by (0, 1, u), where L x 2 dv{x) = 1. Moreover, assume that J R \x\ dv{x) < oo 
and f Ro x 4 dv(x) < oo. Let Y t = j* V2\e- X ^' s) dW s and Z t = ff [0t]xRo V2Xe- x{t - s) xdN(s,x), so, Y t is a 
Wiener O-U process and Z t is a Poisson (pure jump) O-U process. If ht(s) = V2\e- x{t - s) then the 
double O-U process mentioned above is just Y t + Z t = h(h t ). Now, define h[°\s,x) — h t (s)l{ x=0 }(x) and 
h^^SjX) — ht(s)l{ X9 Lf)y(x), then Yt = Ii(h\ ) and Zt = Iiihf^). Notice that due to the normalization of 
the Levy triplet we have that C(t,s) 



LI 



ht ,hs ) ■ The goal of this subsection is to 

' LI 



prove that F T = T~i Y t Z t dt ^ Z ~ 7V(0, E 2 ) as T — > oo. 

Since h[°' and ftf' have disjoint supports (and using the product formula dTTJ)), Y t Z t = I 2 (h[ 0) dih^) , 

and by Fubini Ft = I 2 (T~^ j^. hf ^h^dt). Hence Ft lies in the 2 nd chaos. According to Corollary Q] we 
just need to check conditions (pfe]>. d26j) , ((2Tj) , ([28]) and (|29|) . 



Expectation of the First Derivative's Norm: 



First Malliavin Derivative: 

D Z F T = T^ I h{hf ) )h^\z) + h{h { i l) )h { l ) \z)dt 



— Norm of the First Malliavin Derivative: 

\\D z F T \\l, =T^ / ll 1 {hf ) )h^\z) + h{h^)hf\z)M^ ) W 2 dtds 



=T-> / h{h^)h{hf)) U$\hV>) +h(h?>)h{hV) (h?\h(< dtds 

J[0,T] 2 X 1 L l \ , L; 

=T- 1 I (h(h^)h(h^) +h(h < t 1) )h(h^))C(t 1 s)dtd S 



then, 



\D Z F T \\ L 2 <2T- 



[0,T] 
r 4 



I 1 (h[ 0) )l 1 (h ( ~°'>)C(t,s)dtd i 



J[0,T\* |M Ll 



[0,T] 2 

C(t 1 ,t 2 )C(t 3 ,U)dt 



h(h[ 1) )l 1 (h^)C(t,s)dtd( 



Expectation of the First Derivative's Norm: 

Notice that by the product formula (fTT|) we have that 



,(0) h (0)\ 



K,K) 



Ll 



and 

4 



i=l 
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so, 



E 



1=1 



C{t u t 2 )C{HM) + (h{yh ( t °\h ( t l ) ®h£ ) 



L 2 



and 



E 



n^(^) 

Also notice that 



rtf\S 

C(t, s) = / 2Xe-^ t+s - 2u Uu = e - A l*- s l - e- A (*+ s ) < er A l*-l < 1 
Jo 

min{ti,t2,t3,t4} 



L 2 , 



/■min|ti,r2,i3it4^ /• /• 

/ / a; 4 4A 2 e- A (* 1+t2 + t3 +* 4 - 4tl )^(x)dw < A / x 4 du(a 

Jo JR Jin 



L 2 , 



L 2 , 



C(tl,t3)g(t2,*4) + C l (t 1 ,t4)C'fe,*3) < j 



Putting all this together we get, 



E[\\DF T \\l 2 ] =2T- 



f (E\T\h(h^)] +E\]ll 1 (h£ ) )])c(t 1 ,t 2 )C(t 3 ,U)di 

J[0,T\*\ l l=1 J ItJi \) 



<2(4 + A/ x 4 du(x) ) ( T -1 / C{t,s)dtds 

J[0,T] 2 



e- xlt ~ sl dtds 

[0,T]2 
T /-t 



<2 ^4 + A y x 4 dj/(x)j (T" 1 

=2(4 + \f x i dv{x))(2T- 1 [ f e-W-'^dsdt) 

V ^K J \ JO JO J 



<2 4 + A 



/ x 4 du(x) 

Js.n 



All this proves that 



^[IIDFtII^] 3 = 0(1) as T^oo 



• Expectation of the Cube of the First Derivative's Norm: 

— Cube of the First Malliavin Derivative: 

Since hf ] {z) ■ hi 1] (z) = for all z € R+ x K then, 



ID F T 



Norm of the Cube of the First Malliavin Derivative 
Since x ■ hf\z) = for all z = (t, x) e M+ x M ther 

3 



x\ADF T f) <T"i 



(0)1 



(i) 



i 2 , 
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Expectation of the Cube of the First Derivative's Norm: 
Notice that by Holder we have, 



r 3 



E 



i=l 



<(|C(t 1 ,i 2 )| + s[(/ 2 (/ l w®4°))) 2 ]^£;[(7 1 (<))) 2 ]' 

<(Wi,t 2 )| 



L 2 . 



■,(0) 

*t3 



i 2 



< 2 



and 



u I \ /.min{ii ,i 2 ,t 3 } /■ 

[T \h { t 1} ) =/ {2X)ie- x ^ +t2+t3 ~^du \x\*du(x) 

i=l' / L 2 ^0 7M 



<- 



.2V2A 



(./to 



\x\ dv(x) I e 



~A(ti+t + t 3 — 3 min{ti,t 2 ,t3}) 



Putting all together we get, 

3 / 3 



(i) 



i=l 



\JR / J[0 



< 



e -A(ti+t+t 3 -3 min{ti,t2,t3})^£* 

[0,T] 3 

/ |a;| 3 di/(;r)V~* III e- x{t+s - 2u Uudsdt 
4^(/ Ro |x| 3 ^)) 
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o Jo Jo 
0(T"2) 



All this proves that 



E[(\x\ ,\DF\ 3 j L2 ] -> as T^oo 



Expectation of the Contraction's Norm: 

— Second Malliavin Derivative: 

i-T 



D 2 ZuZ2 F T = T-i I h ( t 1 \z 1 )h ( t °\z 2 ) + h ( t °\z 1 )h ( t 1 \z 2 )dt 
Jo 



Contraction of order 1: 



C(M)<1 



C(t,s)<l 



D 2 F T ®i D 2 F T =T 



1 / ^\ Zl )h^{ Zl )(hf\h^) +$\zx)hV>[z 1 )h$\hV>) dtds 

J[0,T] X ' L l X ' L l 

<T~ X I h^\z 1 )h { p{z l ) + ( Zl )hf \ Zl )dtds 

J\0,T] 
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— Norm of the Contraction: 

-.2 



D Z F T ®i D 2 F T \\ 2 L2 < T- 2 



[ /h {0) h (a) h (0) h (0) \ +(h {1) h {1) h {1) h {1) \ 



dt 



Expectation of the Contraction's Norm: 
Notice that 



/,(!)/,(!) h {1) h {1) 

tl "t 2 >"t 3 "*4 / L2 ~~ \"tl "ts ' *3 "*4 /L2 



so, 



^^ e -A(tl+t2+t3+t4-4min{ti,t 2 ,t3,t4})_ e -A(tl+t2+t3+t4)^ 



(h^hflh^hf}) < 24A / T r r Te- A ' '"'dnhulsdl < 44- 



'[0,T] 4 x " " " ' ^ 

Putting all together we get, 



o Jo Jo Jo 



4T 



E[\\D 2 F T ® 1 D 2 F T \\' L2 ]<2T- 2 I /*.«!,(*) »,(0 fc «) 



[0,T] 4 



;, l "/> 1 " r1t<T~ 2 —--^— 

n t x n t 2 ! re t 3 "*4 / L 2 - ^2 _ 



All this proves that, 



S[||U 2 F T (g)iD 2 F T || i2 ] ->0 as T^oo 



• Expectation of the Squared Second Derivative's Norm: 

— Square of the Second Malliavin Derivative: 



= T" 1 / h[°\z 1 )hf\z 1 )h[ 1 \z 2 )h^\z 2 ) + h[ 1 \z 1 )hi 1 \z 1 )h < i 0) (z 2 )hf\z 2 )dtds 

J[0,T] 2 



Inner Product of the Squared Second Malliavin Derivative: 
Samely as above, x ■ h\°^ — 0, so 



(x,(D 2 F T f) T2 -T- 1 / h[ 0) hW (x^^hW) dtds< [ Ixfdv^T- 1 [ h 

* J[0,T] 2 ^ ' L l JRo J[0,T\ 2 



f ] h^dtds 



— Expectation of the Squared Second Derivative's Norm: 
By the above computations we have, 



E (x,(D 2 F T ) 2 ) L2 
All this proves that, 



< 



E 




2 

L l- 


->• as T -» oo 



Existence of the Variance: 
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Since F T = h(T~^ h^^h^dt) then 



Var[F T 



T 

1 / 



[ hf^h^dt = T- 1 [ ( [ 2Xe'^ 2t ^-^dt) d Sl ds 

J-y- T 2 J[o.t] 2 \J Sl vs 2 J 

=T~ l [ ( e -^l«-«l - e-^ 2T - sl - s ^) 2 d Sl ds 2 

J[0 : T] 2 

=T -l f e ~2A| Sl -s 2 | _ 2e -2A(T- Sl A, 2 ) + e -2A(2T- Sl - S2 ) dsids2 
./fO.Tl 2 



All this proves that, 



Var[Fy] — > i £ (0, oo) exists as T —> oo 
A 



Since all five conditions are met, by Corollary [T] we have that Ft N ~ jV(0, t) as T — > oo. Moreover, 
due to the quantitative property of the inequality, is possible to estimate (from the computations above) that 
the rate of convergence to normality is at least 0(T~3) ; i.e., dw(—?==^,N) = 0(T~i} as T — > oo. This 

rate is similar to the one obtained for the linear functionals of Gaussian-subordinated fields with underlying 
process given by the increments of FBM or the fractional-driven O-U, when H £ (0, |). I believe that the 
right rate of convergence for these processes should be 0(T~^)l 
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